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The  Plastic  Zone  and  Residual  Stress 


Near  a  Notch  and  a  Fatigue  Crack  in  HSLA  Steel* 

by 

W.  H.  Schlosberg0  and  J.  B.  Cohen* 


ABSTRACT 


The  plastic  zone  and  residual  stress  around  a  notch  under  load  and  with 
the  load  removed,  and  around  a  fatigue  crack  (at  the  same  stress 
intensity  factor  as  for  the  notch)  have  been  examined,  with  automated 
X-ray  techniques  and  a  microbeam.  There  is  good  agreement  between  the 
measured  plastic  zone  size  and  Hutchinson's  theory  for  a  work  hardening 
material.  Residual  stresses  exist  well  behind  the  tip,  and  vary  with 
depth,  so  that  measurements  of  crack  closure  on-a  surface  may  not  be 
directly  related  to  closure  s^'ess  (which  samples  the  bulk). 
Instabilities  in  the  dislocation  arrangement  can  be  detected  by 
comparing  T-ray  line  broadening  of  bulk  specimens  under  load,  and  with 
the  load  removed. 
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INTRODUCTION 


It  is  well  known  that  the  stress  singularity  at  a  notch  or  crack 
tip  produces  local  deformation,  and  that  the  associated  plastic  upset 
results  in  residual  stresses  in  and  around  this  region:  The  maximum 
extent  of  the  stresses  and  deformation  is  not  directly  ahead  or  above 
the  crack.  While  there  have  been  many  theoretical  and  expermental 
studies  of  these  phenomena,  no  study  (of  which  we  are  aware)  has 
examined  experimentally  both  the  stresses  and  the  plastic  deformation 
simultaneously,  in  two  dimensions,  around  a  notch  or  crack.  This  is 
the  purpose  of  this  study. 


Rice  and  Rosengren^  and  Hutchinson^. 3)  have  obtained 
theoretical  solutions  for  the  shape  of  the  plastic  zone,  employing  the 
Von  Mises  yield  criterion,  and  allowing  for  work  hardening.  The  former 
authors  developed  the  solution  for  plane  strain,  with  the  shear  stress, 
t,  expressed  in  terms  of  the  yield  shear  stress,  ,  and  the  shear 
strains,  y  and  y^  :  T  „  +  afy/y^)11  . 


Here:  y  "  ,  where  t  and  :  are  stress  and 

strain  components,  and  repeated  subscripts  imply  summation. 

Hutchinson's  solution  is  for  plane  stress,  with  flow  described  with  the 
form: 


In  this  case  the  work  hardening  exponent,  N,  is  the  inverse  of  n  in  Eq. 
1.  Both  solutions  assume  that  the  stress  singularity  near  the  notch  or 
crack  can  be  approximated  by  the  first  term  in  an  asymptotic  series 
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expansion.  Plastic  zones  calculated  from  these  theories  are 
illustrated  in  Fig.  1.  For  plane  strain,  the  maximum  extent  of  the 
zone  moves  closer  to  the  crack  plane  as  work  hardening  increases, 
assuming  a  "butterfly"  shape.  The  zone  size  decreases  with  increasing 
work  hardening  exponent,  for  both  plane  stress  and  plane  strain. 

While  these  calculations  are  applicable  to  monotonic  loading, 
their  application  to  fatigue  requires  some  caution,  because  both 
authors  assume  the  stress  is  proportional  to  strain.  Also,  Rice^  has 
indicated  that  there  may  be  two  plastic  zones  ahead  of  a  fatigue  crack, 
the  outer  one  due  to  tensile  loading,  and  the  inner  one  due  to  reverse 
loading  approximately  one  quarter  the  size  of  the  outer  one. 

Experimental  studies  of  the  plastic  zone  are  summarized  in  Table 
I.  Except  for  the  last  entry, the  delineation  of  the  zone  has  been 
rather  arbitrary,  (and  other  than  the  first  entry)  the  agreement  with 
theory  poor.  Fine  et  al.^D  have  noted  that  agreement  is  good  if  the 
stress  for  zero  hysteresis  in  incremental  strain  controlled  fatigue  is 
employed  in  the  calculation,  rather  than  the  cyclic  yield  stress. 

(This  stress  is  much  less  than  the  cyclic  yield  stress.)  A  number  of 
experimentalists  have  noticed  the  "butterfly"  shape  of  the  zone,  for 
example,  refs.  6,  12,  and  21. 

It  is  well  established  that  the  residual  stress  immediately  ahead 
of  a  fatigue  crack  and  parallel  to  the  aplied  load  is  compressive, 
turns  tensile  at  some  distance,  and  then  oscillates  in  sign,  with 
decreasing  magnitude.  (There  is  no  experimental  study  of  a  notch.) 
Indeed,  it  is  also  well  established  that  the  effect  of  an  overload  is 
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to  increase  this  compressive  stress  in  magnitude  and  extent,  and  to 


decrease  crack  velocity  while  the  crack  is  in  such  a  region/22, 2^ 

Such  stresses  have  been  shown  to  aid  crack  closure/22,2^ 

For  the  most  part,  the  stresses  have  been  measured  only  directly 
ahe.d  of  the  crack,  with  X-rays.  Rice's  calculation^  exhibits  a 
region  of  constant  compressive  stress  immediately  ahead  of  the  crack 
tip  (in  the  reverse  plastic  zone),  but  this  i3  due  to  the  assumption  of 
an  ideally  plastic  material.  Two  finite  element  calculations^’2^ 
indicates  the  stress  decreases  in  magnitude  from  the  crack  tip  with  a 
maximum  value  of  approximately  2/3  oy  .  Some  measurements^ ^7, 22, 27, 28) 
do  indeed  show  the  maximum  stress  at  the  tip,  whereas 
others^®*  17,29,30)  report  the  maximum  compression  ahead  of  the  tip. 
However,  this  might  be  due  to  uncertainty  in  locating  the  tip  with 
respect  to  the  X-ray  beam.  In  all  cases,  the  stresses  are  lower  in 
magnitude  than  predicted  by  Rice.  This  difference  has  been  attributed 
(for  example  in  ref.  27)  to  the  size  of  the  X-ray  beam,  but  this  is 
definitely  not  the  case  in  ref.  17. 

Some  of  these  authors  report  residual  stresses  slightly  behind  the 
crack  tip,  but  usually  attribute  this, again, to  the  size  of  the  X-ray 
beam.  However,  the  beam  was  quite  small  in  ref.  22,  and  furthermore, 
stresses  were  found  well  behind  the  tip;  these  are  probably  a  result  of 
the  dislocations  generated  by  the  propagating  crack. 


The  only  two  dimensional  3tudy  of  residual  stresses  (around  the 
tip  of  a  fatigue  crack)  is  that  of  Allison^0);  however  the  uncertainty 
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in  stress  (  ±70  MPa)  is  quite  high,  only  longitudinal  stresses  were 
measured,  and  the  x^-ray  beam  was  much  larger  than  the  plastic  zone. 


EXPERIMENTAL  PROCEDURES 
Specimens 

The  major  objective  of  this  investigation  is  the  mapping  of  the 
region  around  a  crack  tip.  As  the  changes  in  the  profile  and  position 
of  an  X-ray  peak  could  be  expected  to  occur  in  a  region  less  than  1  mm 
in  size,  a  very  small  X-ray  beam  is  required.  In  a  standard  powder 
diffraction  experiment  with  a  beam  typically  10,000  ^m  x  4,000  p.m,  a 
25  nm  grain  size,  and  with  the  beam’s  penetration  to  this  depth, 

60,000  grains  are  Irradiated.  Special  precautions  are  necessary  to 
assure  reproducibility  of  the  diffraction  profiles,  if  a  beam  the  order 
of  100-200u,m  is  to  be  employed.  Therefore,  a  HSLA  steel  was  chosen 
(Inland  Steel  Co.  No.  328),  because  of  its  inherent  fine  grain  size,  5 
p.m,  so  that  the  beam  samples  5000  grains.  (Tests  on  the 
reproducibility  of  the  diffraction  profile  are  reported  below,  in  the 
section  on  X-ray  measurements.)  Its  composition  is  given  in  Table  II 
and  it  was  obtained  in  the  form  of  3.4  mm  thick  sheets.  Samples  were 
prepared  for:  a)  monotonic  tensile  testing  (to  examine  the  mechanical 
behavior  and  the  effects  of  various  levels  of  plastic  deformation),  b) 
to  determine  the  X-ray  elastic  constants,  c)  with  a  center  notch  for 
fatigue  testing.  Strips  152  x  25  mm  were  sheared  from  these  sheets, 
with  the  long  dimension  parallel  to  the  original  rolling  direction. 

Mill  scale  was  removed  with  a  fly  cutter,  with  a  maximum  cut  of  0.5  mm. 
To  minimize  bending,  specimens  were  given  final  dimensions  with  an  end 
mill  and  a  surface  grinder  (50  urn  cuts  under  a  liquid  spray,  with  13  um 


cuts  in  finishing  passes).  The  geometries  of  the  final  specimens  are 
shown  in  Fig.  2.  When  these  specimens  were  machined,  material  was 
removed  from  one  surface  until  the  intensity  of  the  110  peak  was  a 
fixed  value,  the  remainder  of  the  material  being  removed  from  the 
opposite  face.  This  was  done  to  assure  that  the  preferred  orientation, 
which  was  present  in  all  specimens,  was  the  same  at  the  surface  to  be 
exposed  to  the  X-ray  beam.  Finishing  was  carried  out  with  a  chemical 
etch  (200  parts  30%  HjOg*  !5  parts  48%  HF),  followed  immediately  by  a 
wash  in  methyl  alcohol  and  then  in  water.  This  etching  was  continued 
until  the  110  Kari  -  K&2  doublet  resolution  ceased  to  improve. 

In  experiments  in  which  changes  with  depth  were  studied,  a 
lacquer  was  applied,  except  in  the  area  of  interest,  and  the  above 
chemical  etch  was  employed.  The  thickness  removed  was  measured  with  a 
micrometer,  or  from  a  calibration  of  thickness  vs.  time  in  the  solution. 

The  center  notch  (Fig.  2b)  was  produced  with  an  electric  spark 
discharge,  employing  a  pure  copper  electrode.  Only  specimens  with 
narrow  uniform  notches  within  ±3°  of  the  perpendicular  to  the  tensile 
axis  were  employed.  These  samples  were  lightly  etched  again,  after  the 
notch  was  formed. 

Mechanical  Testing 

Monotonic  stress-strain  curves  were  obtained  on  an  Instron 
machine,  employing  a  "clip-on"  extensometer.  Samples  with  various 
amounts  of  plastic  deformation  were  also  obtained  in  this  manner,  to 
compare  their  X-ray  profiles  to  those  near  the  plastic  zone  of  a  crack 


or  notch. 


To  measure  the  X-ray  elastic  constants  (see  below),  monotonic 
tests  were  conducted  in  the  elastic  range  in  situ,  on  a  diffractometer, 
with  a  small  device  designed  for  this  purpose;^1)  this  is  essentially 
a  set  of  grips  in  a  channel,  which  can  be  separated  by  a  thread  with  a 
fine  pitch,  coupled  to  a  gear  reducer.  One  revolution  corresponds  to 
6.35 ^  m  motion.  A  sample  with  a  mounted  extensometer  was  extended 
first  on  the  Instron  machine  to  obtain  strain  vs.  load,  and  then,  on 
the  diffractometer,  the  extensometer  reading  was  converted  to  stress. 

Fatigue  tests  were  carried  on  a  servo-hydraulic  instrument 
manufactured  by  MTS.  The  (pull-pull)  tests  were  conducted  at  10Hz, 
with  an  R  ratio  of  0.03-0.05.  To  minimize  bending  moments,  the  lower 
grip  was  placed  in  Wood's  metal,  which  was  melted  during  the  mounting 
of  the  specimen.  Crack  extension  was  examined  periodically  with  a  40x 
travelling  microscope.  X-ray  measurements  were  obtained  after  the 
total  length  (crack-plus-notch)  was  approxmately  half  the  width  of  the 
sample. 

In  order  to  compare  the  plastic  zone  sizes  observed  in  this  study 

with  the  various  theories  mentioned  in  the  introduction,  the  yield 

strength  (  a  y),  ultimate  tensile  strength  and  work-hardening 

exponent  (n),  are  needed  for  both  monotonic  and  cyclic  loading.  To 

obtain  n ,  the  continuous  portion  of  the  monotonic  tensile  data  was 

fitted  to  Eq.  1  with  t  ,y  replaced  by c  ,  e  .  The  cyclic  yield  stress 

(32) 

was  defined  from  data  for  another  heat  of  the  same  steel,  scaled 

by  13  pet.  to  compensate  for  a  difference  in  the  monotonic  yield 
stress.  The  results  are  given  in  Table  II. 


X-ray  Measurements 


The  X-ray  source  was  a  Rigaku  rotating  anode  generator  operated  at 
58kv,  lOma,  with  a  fine  focus  filament  (0.1  x  1  mm)  and  a  Cu  anode 
observed  in  point  focus.  With  an  intrinsic  Ge  detector  and  a  single 
channel  analyzer,  Fe  fluorescence  could  easily  be  separated  from  the 
incident  CuF^  wavelength.  A  G.E.  XRD-5  diffractometer  was  modified  to 
fit  this  generator /32)  ancj  a  circular  100  |im  divergence  slit  was 
employed  for  studies  of  the  profile,  with  a  standard  0.05°  receiving 
slit.  This  divergence  was  increased  to  400  p,m  and  the  receiving  slit 
to  0.1°  for  stress  measurements,  because  a  weaker  high-angle  peak  was 
involved. 

At  45°  20  (the  position  of  the  110  reflection)  the  beam  was  250 
p,m  x  100  p.ra,  (sampling  some  5000  grains  as  mentioned  above).  The  peak 
intensity  was  5cp3,  with  a  background  of  0.1  cps,  and  a  (peak)  width  of 
0.3°  20.  A  sample  was  oscillated  ±1,5°  to  increase  the  sampling,  and 
with  this  oscillation  the  peak  intensity  varied  less  than  15  pet.  at 
different  points  on  a  sample. 

To  align  some  particular  point  on  the  specimen  in  the  X-ray  beam, 
the  following  procedure  was  developed.  While  observing  under  a  low 
power  microscope,  a  thin  phosphor  dot  100  u®  diameter  was  applied. 
The  specimen  was  then  placed  on  the  diffractometer  in  a  mount  that 
could  be  displaced  in  x  and  y  directions  parallel  to  the  face  of  the 
specimen,  by  amounts  as  small  as  50M-m.  These  motions  were  employed 
until  the  maximum  brightness  from  the  dot  (due  to  the  X-rays)  was 
obtained.  (A  series  of  divergence  slits,  1  mm  to  100  rim,  were  helpful 
at  this  stage.)  A  low  power  microscope  on  an  adjustable  bed  attached 
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to  the  diffractometer  was  moved  until  its  oross-hair  was  centered  on 
this  dot.  The  sample  could  then  be  moved  to  bring  any  desired  location 
into  the  X-ray  beam,  by  simply  moving  such  a  location  to  the  cross¬ 
hair  of  the  microscope. 

X-ray  intensities  were  accumulated  by  point  counting,  and 
processed  with  a  minicomputer  controlled  diffractometer.  The  software 
was  designed  not  only  for  data  collection,  but  for  on-line  analysis,  as 
will  be  described  below.  This  software  also  included  a  routine  for 
automatically  aligning  the  sample  over  the  diffractometer's  axis. 

(The  sample  displacement  was  determined  that  minimized  the  differences 
in  lattice  parameter  calculated  from  different  peaks.)  Computer 
interfacing  included  a  60  Hz  signal  from  the  rotating  anode.  If  this 
signal  vanished  due  to  an  inadvertent  shut  down  as  a  result  of  an  arc 
in  the  generator,  all  data  were  saved,  and  a  simple  restart  procedure 
allowed  measurements  to  continue  after  the  generator  was  functioning 
again. 

Analysis  of  the  Data 
A)  Profiles 

(34) 

Fourier  analysis  of  peak  shape,  as  modified  by  Dexhez  and 

(35) 

Mittmei jer  was  employed  to  obtain  information  on  microstrains  and 

mosaic  size.  The  entire  process  of  data  collection  and  analysis  was 
carried  out  on-line  with  a  minicomputer  control  system  based  on  a  DEC 
PDP3-E  computer.  To  minimize  the  well  known  effects  of  truncation  in 
such  an  analysis,  four  precautions  were  followed:  1)  at  least  ten 
values  were  obtained  for  the  profile  above  50  pet.  of  the  maximum 
intensity,  2)  fhis  number  of  points  was  never  less  than  15  pet.  of  the 
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total  number  of  points,  3)  the  region  of  the  profile  extended  Con  each 
side  of  the  peak)  at  least  four  times  the  full  width  at  half-maximum 
intensity,  4)  analysis  was  carried  out  about  the  center  of  gravity  of 
a  peak,  to  minimize  the  sine  coefficients. 


All  data  were  corrected  for  the  Lorentz-polarization  factor  for 
step  scanning,  the  variation  of  the  structure  factor,  the  Debye- Waller 
factor,  and  the  dispersion-corrected  scattering  factor.  Analyses  were 
carried  out  on  a  sine  scale.  The  "hook"  effect  (the  decrease  in 
Fourier  coefficients  at  very  low  harmonic  number)  was  minimized 
following  ref.  36,  and  the  Fourier  coefficients  were  corrected  for 
instrumental  broadening  by  Stokes'  procedure/3^  The  standard  for 
this  latter  correction  was  one  of  the  annealed  specimens.  The 
resulting  Fourier  cosine  coefficient,  An,  of  harmonic  number  n,  can  be 
witten  :(35) 


(3) 


Here  is  that  portion  of  the  coefficient  due  to  mosaic  size, 

Deff,  and  a^  is  determined  from  the  range  of  the  peak: 

^a3 

l  m  (sin0  -sin©  .  ). 

\  '  max  min 


The  value  of  na^  s  L  is  the  length  of  a  column  normal  to  the 

diffracting  planes  over  which  the  beam  is  averaging  the  effects.  Mso,(e^) 

is  the  mean-square  strain  averaged  over  such  a  column,  and  d^^  is  the 

spacing  of  the  (hk£)  planes  producing  the  reflection.  The  multiple 

order  procedure  for  separating  Ans  and  <e*>  involves  determining  An  at 

each  n,  for  two  or  more  orders  of  a  reflection,  ie.  vs  ~s  .  Then  the 

d 

average  mosaic  size  normal  to  the  (hk£)  diffracting  planes  (Def>f)  is 
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obtained  from: 


(*) 


Because  of  the  low  intensities  of  higher  order  peaks  in  this 
investigation  with  a  microbean  (less  than  1  cps  for  the  220  reflection, 
for  example)  it  was  decided  to  employ  the  single  peak  analysis 


developed  by  Mignot  and  Hondot. 


(38) 


na. 


For  small  n,  A^  *  1  -tr —  . 

Tikz 


9  l 

Also,  from  ref.  (5),  <c  >  *  — —  •  where  G  is  a  constant. 

Q  QA3 

Substituting  these  relationships  into  Eq.  3,  Mignot  and  Rondot  showed 


that: 


LlS>ktD^t 


(5) 


-  a  +  nfl  +  BaY  •  (6> 

By  algebraic  manipulation  of  Eqn.  (5): 

Deff-  2»3/[-P  +  <P3 -Av>%]  (7a) 

G3  "  {'8  "<S3  ‘  <7b> 

A  least  squares  solution  of  Eqn.  6  (fora  , 3  , y  )  was  obtained 
with  various  combinations  of  low-order  An  (but  excluding  A0).  All 
solutions  involving  the  first  4-10  coefficients  were  obtained.  Those 
with  negative  were  rejected.  The  remaining  solutions  were  ranked  (by 
the  software)  by  considering  that:  1)  a  should  be  unity,  2)  9  should 
be  the  initial  slope  of  An  vs  n,  3)  the  unbiased  residual  should  be  a 

1  0 


minimum. 


Both  the  single  and  multiple-peak  methods  and  all  corrections  have 
been  fully  implemented  in  the  software.  Errors  in  the  resultant  values  of 
the  raicrostrain  and  particle  size  were  obtained  from  the  variances  and 
co-variances  of  the  Fourier  coefficients,  which  depend  on  the  number  of 
counts  collected  across  a  peak.  The  equations  for  these  are  given  in 
the  Appendix.  The  software  was  written  so  that  counting  over  a  peak 
was  repeated  until  the  root  mean  relative  variance  of  the  first  few 
Fourier  coefficients  (which  are  the  ones  that  are  important  in 
determining  Deff  and  G)  was  an  operator-specified  value, : 


Actually,  the  square  root  of  the  sura  of  the  squares  of  for  the 
reference  and  broadened  profiles  was  employed. 

Such  automation  does  more  than  minimize  manual  operations.  It 
also  minimizes  the  time  to  obtain  a  reasonable  precision.  In  the  past, 
this  type  of  analysis  has  been  carried  out  by  obtaining  the  data, 
plotting  ana  smoothing  it,  subtracting  background  by  hand,  and  punching 
cards  for  a  program  for  a" large  computer  to  perform  corrections  and/or 
the  Fourier  analysis.  In  general,  no  error  analysis  is  possible  and 
data  are  usually  obtained  for  times  considerably  longer  than  needed. 
Some  idea  of  the  error  can  be  obtained  by  repeating  measurements  and 
analysis,  but  this  is  rarely  done.  This  older  procedure  is  still 
necessary  for  very  broad,  weak  peaks  with  low  peak-to-background 
ratios,  but  the  new  procedures  described  here  are  applicable  in  most 


situations.  It  is  possible  to  obtain  the  data  and  a  complete  2-peak 
analysis  with  normal  beams  in  3  hours. 

Comparisons  were  made  of  the  single  and  multiple  peak  procedures, 
employing  data  obtained  in  several  past  studies  in  our  group.  In 
general,  if  the  peak-to-background  ratio  is  large,  and  the  mosaic 
size  is  200-500  A  ,  the  single  peak  method  is  viable;  otherwise  it  is 
not.  In  particular,  values  for  (e^>.  are  very  poorly  determined 
outside  this  range,  although  the  particle  size  is  satisfactory. 
Fortunately,  our  studies  fell  within  these  boundaries.  Some  of  the 
comparisons  we  have  made  are  given  in  Table  III. 

B)  Analysis  for  Residual  Stress 

The  222  reflection  was  employed,  which  occurs  at  1 36°  20  with  the 
CuK^  radiation  employed  in  this  research.  The  maximum  intensity  was  2 
cps  (with  a  background  of  0.1  cps).  Because  this  peak  occurs  at  angles 
somewhat  lower  than  those  commonly  employed  for  stress  measurements 
from  steel  with  Co  or  Cr  radiation  (and  which  are  too  low  in  intensity 
for  thi3  study)  the  peak  position  is  more  sensitive  than  usual  to 
sample  position,  often  the  major  source  of  error  in  stress 
measurements.  Therefore,  particular  care  was  taken  to  be  sure  that  the 
surface  of  a  specimen  was  within  25  ^m  of  the  center  of  the 
diffractometer.  This  peak  does  have  an  advantage  though,  in  that  it  is 
unaffected  by  preferred  orientation  and  the  attendant  elastic 
anisotropy  this  can  cause  strong  oscillations  in  d  spacing  vs  sin2i|; 

(where  <|r  is  the  tilt  of  the  specimen  from  the  parafocussing  position); 
it  is  from  the  slope  of  such  a  plot  that  the  stress  is  calculated. 
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Four  v  tilts  were  employed  initially,  from  0°  to  45°,  but  this  was 

reduced  to  three  when  it  was  found  that  ndn  vs  sin2  if  was  indeed  quite 
linear  (correlation  coefficient  >  0.93).  At  each  tilt,  the  same  3° 
oscillation  employed  in  studies  of  the  profile  was  also  used.  A  multi¬ 
point  parabola  was  fit  to  the  top  15  pet.  of  the  peak,  following  ref. 
(44).  The  entire  process  was  automated,  as  described  in  this 
reference.  Statistical  counting  errors  and  geometric  errors  (also  see 
ref.  ^^))  were  evaluated  in  the  software  and  were  typically  a  total  of 
±  20MPa,  which  was  confirmed  by  repeated  measurements. 


To  obtain  the  appropriate  X-ray  elastic  constants,  a  tensile 
specimen  (Fig.  2)  was  mounted  in  the  small  tensile  jig  described  above. 
The  stress  was  kept  below  2/3  g y  ,  to  minimize  plastic  deformation  at 
the  surface.  The  slope  of  "dn  vs.  sin2  t  was  obtained  for  various 
stresses.  Now: 


S2  S2 
where  -f  is  the  effective  elastic  constant.  From  A  vs^  ,  —  was 

obtained,  so  that  Eqn.  9  could  then  be  employed  for  specimens  with 

5 

unknown  stresses.  The  value  of  was  5.08(1.26)  x  10”^  MPa.-1  With 

the  bulk  elastic  constants  in  ref.  46,  an  average  value  of  this 

(47)  -6 

constant  for  constant  strain  and  constant  stress  gave  4.98  x  10 

MPa.-1  The  experimental  value  was  employed  for  all  stresses  reported 
here. 
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RESULTS 


A  Working  Definition  of  the  Plastic  Zone 

The  results  of  the  Fourier  analysis  of  peak  shape  after  tensile 

elongation  are  presented  in  Table  IV.  Distinct  changes  in  mosaic  size 

and  microstrain  occur  at  0.1  pet.  permanent  offset,  after  which,  and  until 

necking  begins,  these  quantities  are  approximately  constant.  It  has 
(48  49) 

been  shown v  *  that  the  mosaic  size  and  microstrain  are  related  to 
the  dislocation  spacing,  and  therefore  to  the  dislocation  density,  p  . 

From  the  mosaic  size  "Qq"  can  be  calculated: 


12  <e»  j>/b2 


where  b  is  the  Burger's  vector. 

The  limit  of  this  analysis  in  this  study,  due  to  instrumental 

* 

broadening  of  the  110  peak,  corresponds  to  a  mosaic  size  of  3500A  . 
Therefore,  Table  IV  implies  a  change  in  dislocation  density  at  0.1  pet. 
offset  from  ~  8  x  1012  m-2  to  1.6  x  1014  ra-2.  Accordingly,  we  have 
chosen  Deff  =  2500A  (9=  1.6  x  lO1^  m“2)  to  delineate  the  plastic 

zone.  This  region  is  shaded  in  several  of  the  subsequent  figures.  It 
is  worth  repeating  that  the  size  of  the  X-ray  beam  on  a  sample  was 
always  considerably  smaller  than  this  zone. 


Residual  Stresses 

One  sample  2  mm  thick,  with  a  5.7  mm  notch  and  a  root  radius  of 
64  p.  m,  was  stressed  at  308  MPa  (—  =  0.52)  so  that  Kt  (calculated 

°y  1 


following  Paris  and  Sih/^O)  including  the  correction  for  the  finite 
width  of  the  sample)  was  30.2  MPa  m1/2.  The  sample  was  unloaded, 
loaded  again,  unloaded  and  loaded  once  more.  It  requires  about  10^ 
cycles  at  this  load  to  initiate  a  crack  from  a  notch  of  this  kind,  so 
this  procedure  produces  a  plastic  zone  ahead  of  the  notch,  without 
propagation.  The  longitudinal  and  transverse  stresses  measuring  under 
load  are  shown  in  Fig.  3.  after  removing  the  load  in  Fig.  4,  and  , 
after  etching  to  remove  one  quarter  of  the  thickness  to  the  center  of 
the  specimen.  Fig.  5.  [This  etching  was  carried  out  unly  in  the 
vicinity  of  the  notch,  by  masking,  as  indicated  in  the  procedures. 
Therefore  no  appreciable  relief  of  stress  due  to  thi3  removal  of 
material  was  expected,  and  no  corrections  were  applied  to  the  data.] 

In  Fig.  3,  it  can  be  seen  that  about  0.8  mm  above  the  notch,  the  measured 
longitudinal  stress  under  load  was  (within  experimental  error)  the 
applied  stress.  Note  also  that  after  removing  the  load  the  compressive 
stress  near  the  notch  tip  is  large  at  the  surface  but  decreases 
considerably  with  depth. 

A  second  specimen,  2.05  mm  thick,  with  a  notch  4.8  mm  long  and  a 
root  radius  of  29  ^  m  was  subject  to  fatigue  at  a  maximum  stress  of  208 
MPa  with  a  stress  range  of  200  MPa,  for  70,000  cycles,  after  which  a 
crack  had  grown  3.5  mm  from  both  ends  of  the  notch.  The  value  of  Kj 
was  17.7  MPa  m1/2  at  the  beginning  of  the  test,  and  31. 3  MPa  m1/2  at 
the  end.  This  latter  value  is  quite  close  to  the  value  employed  above 
for  the  notched  specimen.  The  measured  residual  stresses  are  given  in 
Fig.  6.  The  maximum  value  is  less  than  the  value  of  2/3  <7y  predicted 


Of  particular  interest  is  the  presence  of  residual  stress  behind 
both  the  notch  and  fatigue  crack,  and  the  fact  that  the  stresses  extend 
well  beyond  the  plastic  zone. 


Analysis  of  Peak  Shape 

Typical  errors  in  particle  sizes  (which  were  the  order  of  1200- 
300QA)  and  microstrains  were  25  pet.  The  Figs.  7,8,9  exhibit 
dislocation  densities  calculated  at  various  locations  for  the  sample 
with  a  notch,  and  for  the  fatigued  sample.  The  value  shown  is  the 
square  root  of  the  product  of  Eq.  10  and  11,  that  is,  the  average  of 
the  two  values.  Due  to  the  errors  in  D  and  the  uncertainty  in  these 


values  is  50  pet.  The  density  was  usually  smaller  by  a  factor  of  two 


to  three  when  calculated  from  the  mosaic  size,  implying  that  the 
dislocations  are  clustered.  It  is  particularly  interesting  that  the 
density  immediately  ahead  of  the  notch  (in  the  statically  loaded 
specimen)  increases  when  the  load  is  removed.  This  result  implies  that 


dislocations  move  away  from  tangles  and  walls  on  unloading.  The 


unstable  nature  of  dislocation  arrays  in  the  early  stages  of  fatigue  is 
well  known, Also,  the  density  Just  ahead  of  the  fatigue  crack  is 


higher  than  ahead  of  the  more  blunt  notch,  and  there  is  more  clustering 
of  dislocations  ahead  of  the  fatigue  crack;  the  opposite  is  true  above 
the  crack.  These  patterns  are  in  general  agreement  with  Mugrabi's  TEM 
studies/®1  ^  The  dislocation  densities  just  ahead  of  the  fatigue  crack 


found  in  this  study  are  of  the  same  order  of  magnitude  as  those  found 

by  Yokobori  et  al/1^’1®^  with  the  Hirsch  microbeam  technique  applied 
to  a  low  carbon  steel.  However,  they  reported  a  decrease  of  two  orders 

of  magnitude  in  density  200  n  m  ahead  of  the  crack*  We  see  much  less 
variation. 
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(19) 

Yokobori  and  Sato  examined  the  density  of  dislocations  near  a 

crack  (in  a  low  carbon  steel)  at  various  positions  below  the  surface, 

and  found  little  change  up  to  400  H-ra  from  the  surface.  On  the  other 

(52) 

hand,  Pangborn  et  al.  employing  a  smooth  aluminum  fatigue  specimen, 
reported  a  d ecrease  by  a  factor  of  three  100  p.m  below  the  surface  (followed 
by  an  increase  again  at  greater  depths).  We  made  measurements  of  the 
peak  breadth  at  four  positions  near  the  fatigue  crack  at  50,  130,  190 
and  250  |im  below  the  surface,  and  there  was  no  noticeable  change  in 
broadening.  It  seems  clear  that  near  a  fatigue  crack,  the  dislocation 
density  does  not  vary  appreciably  with  thickness.  This  is  actually  to 
be  expected  from  Pangborn's  results,  which  indicate  that  failure  occurs 
when  the  dislocation  density  in  the  interior  rises  to  that  near  the 
surface.  Such  a  situation  would  be  likely  in  the  plastic  zone  just 
ahead  of  a  crack  or  notch. 

DISCUSSION 

The  specimen  thickness  used  in  this  study  (2  mm),  is  much  less 
than  required  for  plane  strain  conditions  to  dominate;  from  ref.  6, 
this  thickness  would  be  approximately  10  mm  for  the  steel  used  here. 

For  this  reason,  and  also  because  the  X-rays  sample  only  the  near 
surface  regions,  measurements  on  the  face  of  the  samples 
should  resemble  what  is  expected  for  plane  stress.  For  the  sample  with 
a  notch,  the  fact  that  the  residual  stress  pattern  extends  much  further 
normal  to  the  notch  than  ahead  of  it  suggests  that  plastic  upset 
inside  the  specimen,  where  conditions  for  plane  strain  exist,  are 
important  even  near  the  surface. 


The  "butterfly"  shape  ahead  of  the  notch  or  crack  is  clear 
in  our  results  forayy,  but  not  in  the  plastic  zone  itself  (defined 
here  as  equivalent  to  0.1  pet.  plastic  offset  in  tension).  The  regions 
of  residual  stress  need  not  have  suffered  plastic  deformation  but  could 
develop  due  to  upset  in  the  smaller  plastic  zone.  The  stresses  would 
of  course  affect  stress-strain  hysteresis,  and  this  may  be  the  reason 
that  a  very  low  stress  is  required  in  ref.  21  to  calculate  a  "plastic 
zone"  of  the  size  the  authors  measured;  the  zone  delineated  by  stress 
hysteresis  may  actually  be  the  region  of  appreciable  residual  stress. 
With  the  data  obtained  in  this  study  on  yield  stress  and  plastic  zone 
size,  a  direct  comparison  of  calculated  and  measured  plastic  zone  sizes 
is  possible.  This  comparison  is  shown  in  Table  V,  where  it  can  be  seen 
that  the  agreement  is  quite  good  for  the  expected  conditions  of  plane 
stress,  especially  for  the  specimen  with  a  fatigue  crack. 

According  to  Rice^  and  Matsuoka  and  Tanaka^"^  the  plastic  zone 
size  is  ~  1.6  times  the  position  ahead  of  the  crack  where  the  stress 
reverses  sign.  This  value  is  1000  lira,  only  30  pet.  bigger  than  the 
measured  value,  so  this  appears  to  be  a  viable  method  for  estimating 
the  size.  A  reverse  plastic  zone  of  250  would  also  be  expected  in 

this  case.  Although  the  size  of  the  X-ray  beam  employed  in  cur 
experiments  was  almost  half  this  value,  we  could  not  detect  any  unusual 
broadening  or  stresses  very  close  to  the  crack  that  would  suggest  such 
a  region. 

Because  crack  closure  is  affected  by  residual  stresses,  it  is  of 
particular  interest  that  the  stresses  are  much  lower  inside  the 


specimen  than  at  the  surface.  There  are,  of  course,  other  factors  that 

can  lead  to  different  closure  at  the  surfce  and  in  the  bulk,  such  as 

(54)  (55) 

oxidation,  and  the  differences  in  the  stress  state;  all  three 

reasons  complicate  the  relationship  between  crack  closure  measured 

optically  and  closure  stress. 


SUMMARY 

1)  A  quantitative  X-ray  st.tdy  of  the  substructure  and  residual 
stresses  has  been  made  insia?  and  outside  the  plastic  zone  associated 
with  a  notch  and  with  a  fatigue  crack.  A  deformation  corresponding  to 
little  as  0.1  pet.  plastic  offset  in  tension  could  be  detected. 

2)  The  instability  of  dislocation  arrangements  in  such  regions  is 
clearly  indicated  by  changes  in  X-ray  peak  shape  under  load,  vs  load 
removed. 

3)  There  is  good  agreement  between  Hutchinson's  theory  for  plane 
stress  and  experiments  on  the  size  and  shape  of  the  plastic  zone,  for 
both  a  notch  and  a  fatigue  crack. 

4)  There  are  appreciable  residual  stresses  behind  a  notch  or 
crack,  as  well  as  above  and  ahead  of  this  region. 

5)  The  residual  stress  distribution  can  vary  appreciably  with 
depth.  As  a  result  of  this  (and  other  variations  between  the  surface 
and  the  interior)  optical  measurements  of  crack  closure  at  a  surface 
may  not  be  simply  related  to  the  stress  for  crack  closure. 
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APPENDIX 


Errors  in  the  Warren-Averbach  Fourier  Analysis  of  Peak  Shape  Due  to 
Counting  Statistics,  and  Autonation  of  the  Analysis. 

A.  Multiple  Peaks ,  Eq. 

It  is  rare  to  have  more  than  two  orders  of  a  given  reflection 
within  the  observable  29  range,  so  we  deal  with  this  case  here.  This  is 
readily  extended,  if  nore  peaks  are  available.  The  peaks  are  given 
subscripts  1,  2. 


From  Eq.  3,  Aj,  is  linear  vs  1/dg^k£.  Therefore,  by  manipulating 
this  equation  for  each  of  two  peaks  (subscripts  1  and  2),  to  yield  the 
slope  and  intercept  of  An  vs  t/d,  the  variance  (®*  )  can  be  written  as: 

a2 (intercept  -  { (o  (AQl)/di)s  +  (o(Aa2)/d3)3]  (A-la) 


/{l/d3  -  l/d|}3, 

and: 

ea(slope)  -  (o3^!)  +  c3  (A^) )/ [l/d3  -l/d|l. 

The  equation  for  (AQ)  is  discussed  below,  section  C). 

From  ref.  38: 


(A- lb) 


o 3<«£;  ■  (&  <c3)V&  slope)3  a2  (slope) 


+■  (9<«3)^/  8  lntercept)3ca (intercept) ,  (A-2) 

n 


and  hence:  c3(a!)  -  a3 (intercept) , 


o3«c35*)  *  5-<«3>{o3  (slope) /slope3 

•**  g2  (intercept)/ intercept3}. 


(A- 3a) 


(A- 3b) 


The  error  for  the  particle  size  (see  Eqn.  3)  can  then  be  obtained  from 

s 

the  error  in  the  slope  of  a  (weighted  linear  least-squares)  fit  to  Aq 
vs  n,  for  small n  . 
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At  least  four  Fourier  coefficients  of  low  order  (low  n  or  L)  are 
employed  to  solve  this  equation  by  least-squares,  as  for  example,  in 


Ch.  4  of  ref.  39.  However,  the  first  coefficient,  A0,  cannot  be 
employed;  its  value  is  unity  by  definition,  and  therefore  its  variance 
and  co-variance  with  other  coefficients  is  null.  Again,  following 
ref.  38,  and  keep  in  mind  that  the  least  squares  analysis  of  Eq.  6 
yields  ®s  (  »  ),  9a  (y): 

*a(Q)  -  (0)  9a<0>  +  2  If  5^  COV(8,y)  +  *a(Y), 

Here  "COV"  means  covariance. 

we  define  Q  a  -  0*  (^“^Y^for  the  calculation  of  o  (Deff), 
and  Q  a  -  0  -(^-^Y^for  ff(G3).  since  from  Eqs.  5,6: 

D#ff  *  2a3/{-8  +  ©3  -  4y)^}, 

G3  -  daC-B-fe3-^)^/^  a3. 

Therefore,  with  these  definitions  and  Root  *  (  02  -  4y  ^  it  can  be 
shown  that: 

9<Dtff)  -  t(“l  +  0/ROOT)3  o3  (8)  +  4(-l/ROOT)3  ea(Y) 
+4(-l-t*B/R0OT)(-l/R0OT>C0V(B,Y)l^Deff/2a3 

and: 

©(C8)  -  {(-1  -8/R00T)3c3(8)  +  4(1/R0OT)V(y) 

+  4(-l  -0/ROOT) (1/ROOT)COV(0,y)}^ 

*  ds/4tr  a3. 
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(A-4) 


(A-5a) 

(A-5b) 


(A-6a) 


(A-6b) 


The  least  squares  analysis  fora  ,y  in  Eq.  6  require  a  knowledge  of  the 
variance  of  Fourier  coefficients  As  (see  ref.  16)  as  does  the  error 
analysis  for  two  or  more  peaks  in  part  A.  Therefore  we  turn  now  to 
this  variance. 

C.  Variance  of  the  Fourier  Coefficients 

Vilson<40‘42)  has  derived  equations  for  the  variances  of  the 

m  3 

Fourier  cosine  coefficients  is  (Aq)  )  and  sine  coefficients  (®  (B  q)  )  of 
a  Bragg  peak,  as  well  as  their  co-variances.  Be  ignored  certain  small 
terms,  which  we  include  here,  as  the  calculations  are  to  be  carried  out 
on  a  computer.  For  details  of  the  derivations,  the  reader  is  referred 
to  Wilson's  papers  and  ref.  32. 


(A-9) 


of  A  may  be  expressed  as: 


a  KY 

»  (An)  -  drr) 

*X V 


e3dk). 


(A- 10) 


(40) 

For  fixed-time  measurements: 

»a(2k>  -  yr. 

where  T  is  the  measurement  time  per  point. 

From  Eq.  A-9  and  A-10,  with  L0  the  background  corrected  integrated 

intensity:  g»  tok  +  -  g]  cos(4nnk/R) 

o 

+  gCcos(4rrak/R)  -  2AQL[l^  -  g]  coa(2rmk/R) 

-  2AagZcos(2rmWR)  +  .  (A-ll) 

The  individual  terms  in  the  previous  expression  may  be  rewritten  as: 
first:  (Lc  ♦  Rg)/2, 

second:  Le  A^/2, 

fourth:  -2^1^  • 

sixth:  Ag(L0  *  RS>*  (A-12) 

Combining  these  terms  allows  equation  (A-11)  to  be  rewritten: 


°3(An>  *  tF  (Au<Rg  "  Lo)  +  U2nLo]/2  +  [Lo  +  ^/2 
o 


+  gC^Scos (4rmk/R) - 2A  Lcos (2nnk/ R) ] 1  . 


(A- 13) 


A  similar  analysis  starting  with  the  definition  of  the  covariance 
between  the  nth  and  ath  Fourier  coefficients: 


C0V(A 


/&*  *a\ 


<V 


(A- 14a) 


gives  the  following  result: 


C0V<VAJ  "  tF  ^Vm(R*  *  Lo)  +  CAo+«Lo]/2 
o 


IB  \  a 

a\  _a , 


°a<v-i  *(1>- 


(A- 15b) 


giving: 


•a<V  -  HF  -  V  "  U^/z+CV**]/* 


*C%gE«on(Annk/R) 

+2B«*®k  ”  CL>/*23t*i||(2T,n,l^Il)31  •  (A-15e) 

Equation  (A— 12)  can  ba  used  as  tha  criterion  for  determining  the 
tiaa  for  data  collection,  as  wall  as  for  the  analysis  of  errors.  Given 
that  the  initial  tiaa  for  data  collection  which  is  specified  by  the 
user,  T,  is  long  enough  to  measure  a  statistically  significant  number 
of  counts,  the  total  tlae  of  data  collection  can  ba  predicted  froa  the 


expression: 

•»*  T  (A-16) 

(A>  " 


where  c®  (A)  is  the  variance  calculated  for  a  measurement  for  time  * 
a 

and  o  p  (A)  is  the  desired  variance;  T  is  the  required  counting  time. 
In  practice  Che  Fourier  coefficients  will  vary  slightly  as  a  function 
of  the  counting  time  and,  therefore,  the  predicted  total  time  may  prove 
insufficient.  Since  the  process  is  iterative,  the  sequence  of  steps 
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is: 

1)  measure  quickly, 

2)  calculate  Fourier  coefficients  and  variances, 

3)  calculate  additional  counting  time. 

This  sequence  may  have  to  be  repeated  a  third  or  even  a  fourth  time. 

The  method  of  specifying  the  desired  variance  should  take  into 
account  that  a  number  of  Fourier  coefficients  are  needed  to  determine 
the  particle  size  and  root-mean-square  strain.  One  way  to  accomplish 
this  is  to  define  the  allowable  error  as  the  root-mean  relative 
variance: 

«d<a) 

However,  since  only  the  initial  Fourier  coefficients  are  used  to 
determine  the  strain  and  particle  size,  '  can  be  limited  to  Include 
coefficients  up  to  an  arbitrary  maximum.  In  our  case  coefficients  up 
to  a  column  length,  L  =  naj,  of  200  A  with  an  arbitrary  maximum  of  n  *  5 
were  used.  Since  Aq  is  unity  by  definition,  it  in  not  included  in  this 
calculation.  If  Fourier  coefficients  are  known  for  both  the  reference 
and  the  broadened  profile,  the  square  root  of  the  sum  of  the  two  root- 
mean  relative  variances  is  a  suitable  estimate  of  the  root-mean 
relative  variance  of  the  Stokes  corrected  profile.  This  was  used  in 
this  study  to  determine  the  counting  times. 

D.  Features  of  the  Program 

1.  An  initial  dialogue  with  the  operator  requests  pertinent 
information,  such  as  the  appropriate  equation  for  the  Lorentz- 
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polarization  factor,  absorption  coefficients,  scattering  factors, 
oscillation  range,  29  limits,  deadtime,  wavelength,  preset  time  or 
count,  29  range  of  peak  and  step  interval  (which  may  be  different  in 
different  regions  of  a  peak),  and  percent  error  desired  in  the  Fourier 
coefficients.  Input  information  which  varies  with  29  is  fit  with  a  cubic 
spline  function. 

2.  As  a  peak  is  analyzed,  various  facets  of  the  analysis  are 
printed  and  plotted,  to  allow  the  operator  to  change  items,  or,  for 
example,  to  choose  a  different  set  of  Aq  in  the  one  peak  analysis.  The 
output  Includes  the  particle  size  and  strain  and  the  associated  errors. 

For  further  details,  see  ref.  32.  A  program  listing  as  well  as  a 
user's  manual  are  available  from  the  second  author. 
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TABLE  V 

Comparison  Of  Measured  Plastic  Zone  Sizes  with  Those 
Calculated  From  Theoretical  Models. 


Measured 


ANGLE  HELATIYI 

TO  CRACE  PLANS, 

0  (•) 

• 

*  (S)  -  ^ 

r(d) /Ilj/Gj  ) 

r(S) 

jun) 

Statically  loaded 

0 

0.19 

500 

sample*!  using  ) 

90 

e.30 

800 

Fatigued  sample 

0 

0.25 

700 

(using  j&tj  } 

90 

0.18 

500 

Fatigued  sample 

0 

0.23 

700 

(using  I  ) 

X  max 

90 

0.15 

500 

dork  hardening 

0 

0.33 

930 

600 

plane  stress 

90 

0*21' 

Work  hardening 

0 

0.00 

0 

plane  strain 

90 

1 

j 

0.2* 

680 

(*)  (I./Cj)8  *  2.82  mm  used  tor  calculated  values 


FIGURE  CAPTIONS 


Fig.  1  Plastic  zone  boundary  for  a  work  hardening  material  subjected  to  mode 
I  loading.  Poisson's  ratio  ■»  0.3.  Zone  for  plane  strain  calculated 
from  ref.  1,  that  for  plane  stress  from  ref.  2,  3. 


Fig.  2 
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Fig.  3  Static  load  applied,  sample  with  notch.  o/oy“  0.52,  Kf  30.2  MPa*  nr*. 

Plastic  zone  shown  shaded  $diere  Dgff  >  2500A) :  a)  longitudinal  stress 
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Fig.  4 


Longitudinal  stress,  measured  after  static  load  removed  for  sample 

with  notch  o/Cy*  0.52,  ■  30.2  MPa* nr*.  Plastic  zone  shown  shaded. 


Fig.  5  After  etching  one  quarter  of  the  way  to  the  center  of  sample  with 

notch,  load  removed,  o/Cym  0.52,  30.2  MPa*nT*.  Plastic  zone 

shown  shaded,  a)  g  ,  b)  o 

yy  xx 

Fig.  6  After  70,000  cycles  in  pull-pull  fatigue,  load  removed.  Ao/ff'  "  0.45, 

A K_*  31.3  MPa'tn3.  Plastic  zone  shown  shaded,  a)  «  ,  b)  c  • 

i  yy  xx 


Fig.  7 


Dislocation  density  (X 
applied.  o/Oy  ■  0.52, 


10 

*1 


^m“2)  for  sample  with  notch,  static  load 
*  30.2  MPa*  nr*.  Plastic  zone  shown  shaded. 


Fig.  8 

Fig.  9 


Dislocation  density  (X  10 
removed.  a/Cy  *  0.52, 


”^m"2)  for  sample  with  notch  static  load 
*  30.2  MPa*oT*.  Plastic  zone  shown  shaded. 


•  13 

Dislocation  density  (X  10  m"2)  for  sample  with  fatigue  crack  after 
70,000  cycles  in  pull-pull  fatigue,  load  removed,  fa/tiy  *  0.45,  AKj  * 
31.3  MPa*m^.  plastic  zone  shown  shaded. 
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The  plastic  zone  and  residu. .1  stress  around  a  notch  under  load  and  with  the  load 
removed,  and  around  a  fatigue  crack  (at  the  same  stress  intensity  factor  as  for  the 
notch)  have  been  examined,  with  automated  X-ray  techniques  and  a  microbeam.  There 
is  good  agreement  between  the  measured  plastic  zone  size  and  Hutchinson's  theory  for 
a  work  hardening  material.  Residual  stresses  exist  well  behind  the  tip,  and  vary 
with  depth,  so  that  measurements  of  crack  closure  on  a  surface  may  not  be  directly 
related  to  closure  stress  (which  samples  the  bulk).  Instabilities  in  the  dislocation 
arrangement  can  be  detected  by  comparing  X-ray  line  broadening  of  bulk  specimens 
under  load,  and  with  the  load  removed. 


",,1473 


(PAGE  I ) 


S/N  0101. 807.6801 


Secuntv  Classification 


1 

mi 


Unclassified 


Security  Classification 


KEY  AOROS 


stresses,  residual  stresses,  fatigue,  plastic 
closure. 


DD  ,”',“..1473  <=*«> 


1  NOV  •• 

(PAGE  2) 


